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ABSTRACT 



Context. The late A and F-type y Doradus (y Dor) stars pulsate with high-order gravity modes (g-modes). The existence of different 
evolutionary phases crossing the y Dor instability strip raises the question whether pre-main sequence (PMS) y Dor stars exist. 
Aims. We intend to study the differences between the asteroseismic behaviour of PMS and main sequence (MS) y Dor pulsators as 
predicted by the current theory of stellar evolution and stability. 

Methods. We explore the adiabatic and non-adiabatic properties of high-order g-modes in a grid of PMS and MS models covering the 
mass range 1.2Mo < M, < 2.5Mo. 

Results. We have derived the theoretical instability strip (IS) for the PMS y Dor pulsators. This IS covers the same effective temper- 
ature range as the MS y Dor one. Nevertheless, the frequency domain of unstable modes in PMS models with a fully radiative core is 
greater than in MS models, even if they present the same number of unstable modes. Moreover, the differences between MS and PMS 
internal structures are reflected in the average values of the period spacing, as well as in the dependence of the period spacing on the 
radial order of the modes, opening the window to determination of the evolutionary phase of y Dor stars from their pulsation spectra. 

Key words, asteroseismology - stars: oscillations - stars: variables: general - stars: pre-main sequence 



1. Introduction 

' The y Dor stars are variable late A and F-type star s. Their vari- 
, ability was identified as caused by pulsation by iBalona et alJ 
■ d 19941) . and^ave et al. ( 1999) classified them as a new class of 
pulsators and defined the features of this group. These stars are 
pulsating with high-order g-modes in a range of periods between 
approximately 0.3 and 3 days. 

The observational y Dor IS covers a part of the Hertzsprung- 
Russel Diagram (HRD) between 7200 - 7700 K on the zero- 
age main sequ ence (ZAMS) and 6900 - 7500 K above it 
! (lHandlei]fT999h . between the solar-hke stars domain and the 
6 Scuti (6 Set) IS. We note that they are located between stars 
with a deep convective envelope (CE) and stars with a ra- 
diative envelope, in the region of the HR diagram where the 
depth of the CE changes rapidly with the effective tempera- 
ture of the star Nowa days 66 stars have been confirmed as 
bona fide y Doradus (Henrv et al. '2007), and thanks to the 
space missions CoRoT (Baglin & Fridlund 2006) and Kepler 
fcilliland et al.i[2010l) , the number of y Dor candidates is rapidly 
increa sing (see e.g . lUvtterhoeven et"al] 20081 iMathias et al.l 
120091 iHareter et al.l lloiO, Balona et al.l 1201 ll) . The observa- 
tional limits of the y Doradus IS have been established last by 
[Handler & Shobbrook (2002) (HS02 hereafter), and in the rest 
of the paper these limits will be adopted to define the observa- 
tional y Dor IS. 




Fig. 1. PMS and MS (respectively full & dashed thin lines) evo- 
lutionary tracks that cross the observational y Dor IS (thick lines 
- HS02). Grey points represent the bona fide y Dor stars from 
iHenrv et al. (2007) . 



iGuzik et"al] (l2000h use the frozen convection approximation 
to propose the modulation of the radiative flux at the base of 
the CE as the excit ation m e chani sm. This mechanism has been 
confirmed by Dup ret et al.l (|2005 | ) using a time-dependent con - 



vection (TDC) treatment (lGabrielll996ifGrigahcene et al.l2005h . 
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Because the depth of the CE plays a major role in the driving 
mechanism of j Dor pulsations, the theoretical predictions of 
stability are very sensitive to the parameter a defining the mean 
free path of a convective element (A = a x Hp, where Hp is the 
pressure scale height ) in the classical mix ing-length treatment 
of cor ivection (MLT lBohm-Vitensel[T958l) . Using TDC treat- 
ment, [Dupret et al.' ('2004') obtained good agreement between 
theoretical and observational j Dor IS for models computed with 
a = 2.00. 

These theoretical works on y Dor stars systematically stud- 
ied MS models. However, as shown in Fig. [1] MS and PM^] 
evolutionary tracks cross the observed y Dor IS. Although the 
time spent by a star of 1.8 M© to cross the IS during its PMS 
evolution is ten times less than the time spent during the MS 
phase in the same region of HR diagra m, recent photo metric ob- 
servations of young clusters (NGC 884, S aesen et al. '2()10) have 
revealed the presence of y Dor candidates that, given the age 
of the cluster, should be in the PMS phase. A str ong effort has 
also b een made to find PMS pulsators in NGC2264 (IZwintz et al.l 
l2009h . another young open cluster. Moreover, the PMS/MS sta- 
tus of H R 8799, a y Dor variable hosting four p lanets (or brown 
dwarfs) (iMarois et al.l[2008 iMarois et al.ll20ld) , is sti ll a matter 
of debate (iMova et alJl2010l]Moro-Ma rtm et al] |2010h . 
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Fig. 2. Evolutionary tracks for different models crossing the 
y Dor observational IS (thick grey lines - HS02). Left panel: 
1 .8Mo evolutionary track whose PMS (full line) and MS (dashed 
line) phases intersect inside the IS (circle: PMS model - triangle: 
MS model). Right panel: evolutionary tracks for 1.9 and 2.1Mq 
showing the same HR location for a model of I.IMq in the PMS 
phase (circle) and a model of 1 .9Mq at the end of the MS phase 
(triangle). 



It is then timely to theoretically study the seismic properties 
of PMS y Dor in order to derive possible differences between 
their spectra and those of y Dor in the MS phase. This is the aim 
of this paper, which is structured as follows. In Sect. |3]we com- 
pare the internal structure of MS and PMS models computed as 
described in Sect.|2] The effects of these structure differences on 
the properties of the adiabatic frequency spectra of PMS and MS 
y Dor are analysed in Sect. |4] In Sect. |5] we focus on the exci- 
tation and damping of y Dor pulsations in PMS models and on 
the differences between PMS and MS non-adiabatic quantities. 
A summary is finally given in Sect.|6] 



2. Stellar models 



Stellar models were computed with the stellar evolutionary 
code CLES (Scuflaire et al. 2008b) for masses between 1.2 and 
2.5Mq, initial helium mass fraction Yq = 0.28, a nd initial metal 
mass fraction Zq = 0.02. We adopted the AGS05 (lAspl und et al.l 
12005!) metal mixture, a nd the correspondin g opacity tables were 
computed with OF|j ( Badnell et al.l l20()5h facil ities and com- 
pleted at low temperature (logT < 4.1) with iFerguson et aP 
(l2005l) opacity tables. We used the OPAL2001 equation of 
state (^Ro gers & Navfonovll2002l) and the nucle ar reaction rates 
from NACRE compilation ( Angulo et al.lll999h . except for the 
''*N(/7, y)'^0 nuclear reaction, for which we adopted the cross 
section from Formicola et al. (2004). Surface boundary condi- 
tions at T - T ^s were provided by ATLAS model atmospheres 
dKurucz" 19981. Convection was treated using the MLT formal- 
ism (Bohm-V ite"nsell958h with an MLT parameter a - 2.00. We 
considered models without and with convective core overshoot- 
ing, with an overshooting parameter Cqv = 0.20 and the distance 
of extra-mixing given by do^ - ^ov x inin(rcc. Hp), where r„. is 
the convective core radius. 



3. Internal structure of PMS and MS models 

During the Hayashi track evolution, the global contraction of the 
star leads to the increase in density and temperature. This tends 
to decrease the opacity and a radiative core develops and grows 
in mass as the star evolves. When the central temperature is high 
enough (~ 1.7 10^ K), the nuclear CN subcycle starts. As the 
'^C(/:', y)'-'N(yS^v)'-'C(/?, 7)''*N nuclear reaction is very sensitive 
to the temperature (oc r'^), a convective core (CC) appears in 
the star. The mass fraction of this CC changes as the star evolves 
toward the ZAMS and for stars slightly more massive than the 
Sun, this CC remains during the MS. 

We see in Fig. [1] that the evolutionary phase at which a PMS 
star crosses the y Dor IS depends on its mass. While lower 
mass stars have already developed a CC, the more massive stars 
still have a contracting, chemically homogeneous radiative core 
when they cross the IS. We recall that the PMS phase is a brief 
evolutionary phase compared to the MS one, especially if we 
only consider the PMS phase before the onset of nuclear reac- 
tions. 

The properties of high-order g-mode spectrum are deter- 
mined by the matter stratification in the star, which is described 
by the Brunt-Vaisala frequency 



1 dlnP dlnp 



dr 



r, dr 



(1) 



with g the local gravity, p the local density, P the local pressure, 
r the local radius, Fi the first adiabatic exponent, Vad and V are 
respectively the adiabatic and stellar temperature gradients, 
the mean molecular weight gradient, and ip and 5 are the par- 
tial derivatives of density with respect to yU and temperature, re- 
spectively. It is possible to highlight the direct link between the 
g-mode spectrum and the internal structure of the star from the 
first-order asymptotic theory (Tassoul 1980). In this appro xima- 



' We consider as PMS models those before the onset of the central 
H-buming at equilibrium. 
- http://opacities.osc.edu/ 
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tion, the period of a high-order g-mode of degree ^ in a star with 
a CC and a CE is written as 



V?(7TT)X:' ^dr 



{2k + 1) , 



(2) 



with k the mode radial order, ro and r\ the limits of the ^-modes 
cavity defined by ci? <k ,S'\, where w is the mode frequency 
and Si the Lamb frequency for ^-degree modes. 

Eq. ^ makes the role played by in the central regions 
evident in the determination of j Dor oscillation properties. 
Because of the differences between PMS and MS stellar struc- 
tures, we expect differences between the seismic properties of 
stars in these evolutionary phases. 

Fig-Epresents the evolutionary tracks of LSM© (left panel), 
1.9, and 2.IM0 (right panel) models. PMS and MS 1.8 Mq evo- 
lutionary tracks cross each other at different points of the HRD, 
in particular inside the y Dor IS. The comparison of two mod- 
els at the same location in the HRD allows us to eliminate the 
effects of different effective temperatures and radii on the stellar 
structure. 

In Fig.[3](top panel), we plot the profiles of the two 1 .8M0 
models. Because of the same mean density and T^s, PMS and 
MS Brunt-Vaisala frequency profiles show similar behaviour in 
the outer layers, and the bottom of the CE is located at the same 
depth in these two models. Both profiles also present a bump 
in the inner layers, which is due to the density distribution. In 
the PMS model, the onset of the i2c(/7, y)i3N(/3+v)i^^C(/:', 7)^^^ 
nuclear reaction has already led to the development of a small 
CC. The MS model, in addition to a larger CC, presents a sharp 
variation in the Brunt-Vaisala frequency profile due to the mean 
molecular weight gradient let by the CC during the MS evolu- 
tion. The profile of a PMS star is smooth, even if a small CC 
already exists. 

As above mentioned, more massive PMS models crossing 
the 7 Dor IS are in a less evolved phase than lower mass ones. In 
order to have MS models with low enough effective temperature 
and high enough luminosity to meet more massive PMS tracks, 
we have to consider models with CC overshooting (ffov = 0.20). 
The intersection between 2.1Mq and 1.9Mq evolutionary tracks 
in the y Dor IS occurs when the PMS 2.1Mq model has a ra- 
diative and quasi chemically homogeneous structure while the 
1.9Mq one is at the end of its MS (Fig. |2]- right). Compared 
with the 1.8 Mq models described in the previous paragraph, the 
2.1 Mq PMS model is less evolved and still has not developed a 
CC; on the other hand, the 1.9 M© MS model is more evolved 
and presents a larger V^. Moreover, although these models have 
the same radius, they do not have the same gravity, so the den- 
sity profiles in the outer layers and the depth of the convective 
envelope are different (Fig.|3]- bottom panel) 



4. Adiabatic analysis - Period spacing 

PMS and MS stractures in the y Dor region differ particularly in 
their central layers. We performed an adiabatic oscillation study 
of our grid of stellar models in order to link these differences 
to asteroseismic quantities. The adiab atic oscillation frequen cies 
were computed with the code LOSC (iScuflaire et al.ll2008ah . 



o 




o.a 0.4 0.6 0.1 
r/R. 
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Fig. 3. Propagation diagram. PMS (black), MS (grey) A^ (full 
lines), and 5^=1 (dotted lines) frequencies for models with the 
same mass (top panel) and for models with different masses (bot- 
tom panel) presented in Fig.|2] The thick horizontal lines repre- 
sent the propagation zone for typical y Dor g-modes. 



4.1. Evolution of the period spacing in the first order 
asymptotic approximation 

From Eq. (|2]l we obtain the expression for the asymptotic period 
spacing between two g-modes with consecutive radial orders and 
the same ^-degree 



AP = P, 



2jt^ 



k+\ 



(3) 



This asymptotic approximation predicts a constant value of the 
period spacing, independent of the radial order k. 

Fig. ID (bottom panel) shows the evolution of the asymptotic 
period spacing of a 1 .8M0 star from the PMS to the end of the 
main sequence. The variation in the period spacing value during 
the stellar evolution seems strongly linked to the evolution of 
the CC (middle panel). The largest difference between PMS and 
MS period spacing at a given effective temperature is about 0.02 
days and corresponds to the maximum mass fraction of the CC 
during the PMS phase (Fig.|4]- point B compared to point D). 

In Fig. |5] we present the ^ ratio for models A to E (with 
X the normalized radius). For high-order g-modes in 7 Dor pul- 
sators, the limits of the propagation cavity almost coincide with 
the boundaries of convective regions, and the contribution of ^ 
to the AP value comes from the whole radiative zone. Because 
of its large CC, model B has a low J ^t/x value, leading to a 
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Fig. 4. Top panel: Evolution of a I.SMq star in the HRD (PMS 
phase in full black and MS phase in dashed grey) crossing the 
J Dor IS (thick black lines - HS02). Middle panel: Variation in its 
CC mass from the PMS phase to the early MS phase (X^ = 0.69). 
Bottom panel: Evolution of the t - \ modes period spacing as a 
function of the effective temperature of the star from the PMS to 
the end of MS (point E). 



high value of Af (see Fig.|4]i. Even if models A, C, and D have 
very different structures, hence different profiles, the values of 
the integral over the propagation regions lead to similar values of 
Af . The high contribution to J ^c/x from the region in model 
D is within an order of magnitude of the contribution of the most 
central layers of model A. On the other hand, the higher value 
of ^ in region is offset by the lower value in the outer layers 
when model D is compared to model C. 



o 

X 




Fig. 5. ^ ratio as a function of the normalized radius for models 
A, B, C (respectively full, dashed, dotted black lines), D and E 
(respectively full and dashed grey lines). 



Fig.|6]illustrates the evolution of the t - \ period spacing for 
models between 1.5 and 2.3M0. Although the value of period 
spacing for models in the PMS phase is generally higher than 
for models in the MS phase, they approach close to the MS. As 
a consequence, we conclude that the average value of Af cannot 
be used alone to distinguish PMS j Dor from MS ones. 



0.05 

CO 

<] 0.04 




9000 8000 7000 

log 

Fig. 6. Evolution of the i - \ modes period spacing as a function 
of the effective temperature of PMS and MS models for evolu- 
tionary tracks between 1 .SMq and 2.3M0 in the T^s range of the 
y Doris. 



4.2. Oscillatory signature in the period spacing 

iMiglio et al.l ( l2008l) investigated the properties of high-order g- 
modes in MS stellar models and show that the sharp variation 
in the profile at the limit of the CC in an MS model lets 
clear asteroseismic signature: the oscillation of the period spac- 
ing around its mean value. They define SPi^ as the difference be- 
tween the periods of a star showing such a sharp variation in A^ 
and the periods of an otherwise fictitious smooth model with the 
same value of —dr. Assuming the ICowlingI (fT941l) . JWKB 



ff' 

(see e.g. lGoughll2007h and asymptotic approximations, they de- 



4 



M.-P. Bouabid et al.: Theoretical seismic properties of PMS y Dor pulsators 



rive the following expression of 5Pk for a profile modelled by 
a step function 



^p.ccHilll-^ cos (2^ + 



(4) 



where <I) is a phase constant, v - {jf^ with and respec- 
tively the values of the Brunt-Vaisala frequency at the outer and 
inner borders of the //-gradient region. The local buoyancy radius 
is defined as 



n:' = -dx 



(5) 



Then the total buoyancy radius is 
N 



-dx, (6) 
and the buoyancy radius at the sharp variation is 

'-dx (7) 



with xi the normalized radius at the top of the propagation cav- 
ity, and xo and x^ the normalized radii respectively at the bound- 
ary of the convective core and at the location of the sharp varia- 
tion of the profile produced by . 

Eq. (HI describes the signature of the abrupt change in the 
profile as a sinusoidal component in the period spacing variation 
with an amplitude proportional to the sharpness of the change in 
and a periodicity in terms of the radial order k given by 



Ak 



no' 



(8) 



Fig- 13 (top panel) presents the { - I period spacing as a 
function of the radial order k for the 1 .8M0 PMS and MS mod- 
els described in Sect.|3] The bottom panel presents the A^ profiles 
of these two models versus the relative buoyancy radius. In the 
case of the PMS model, since is not significant, no sharp 
variation in A^ appears in the radiative region, making AP almost 
constant. On the contrary, the sharp feature in the Brunt-Vaisala 
frequency profile of the MS model introduces an oscillatory vari- 
ation in the period spacing values. For the MS model, we see 

in Fig. I2] (bottom panel) that value is about 5, which is 

the mean number of radial orders between two consecutive min- 
ima of the period spacing value (Fig]?] - top panel). This oscil- 
lation does not show a sinusoidal behaviour as expected from 
Eq. because the approximations used to derive this expres- 
sion are only valid for small variations relative to the smooth 
model. Nevertheless, we verify the relation between the oscil- 
lation periodicity in terms of k and the buoyancy radius of the 
sharp variation (11^). 

Therefore, the deviation of the period-spacing behaviour 
from the constant value expected from the first-order asymptotic 
approximation can be used to distinguish between PMS and MS 
7 Dor pulsators. In fact, the periodicity of the variation of AP 
is directly linked to the location of the chemical gradient at the 
border of the convective core and the amplitude of this variation, 
and its dependenc e on k give information on the A^ profile (see 
iMigUo et alJl2008l for details). 




-10 -15 




Fig. 7. Top panel: Period spacing structure for { - I modes as 
a function of the radial order k of the modes for PMS (black 
triangles) and MS (grey squares) models with the same mass 
(M = I.8M0). Bottom panel: The Brunt-Vaisala frequency ver- 
sus no/flv of these I.SMq PMS (black) and MS (grey) models. 
The dashed vertical line indicates the location of 2ii for the MS 
model. 



5. Non-adiabatic analysis 

This section presents the results of our non-adiabatic oscil- 
lation study. We computed t - I and ( - 2 modes for 
our grid of stellar models using the non-adiabati c code MAD 
with the time-dependent treatment of convection (lDupretll200ll 
iGrigahcene et al.ll2005h . 

In Fig. m we present the theoretical IS for PMS and MS 
^-mode pulsators in the y Doradus region of the HRD. PMS 
high-order g-modes pulsators are expected in this region, and 
the unstable PMS models cover the same effective temperature 
range as the MS y Dor IS. The only difference between both IS 
lies in the extension of the PMS one toward higher luminosi- 
ties. With our set of physical parameters (see Sect. |2}, we cannot 
obtain MS models in the upper region of the HRD where we 
have young massive PMS pulsators, even with CC overshooting 
(ffov = 0.20). 

Fig.|9]shows the £ - \ instability domain versus the effective 
temperature for all PMS and MS models with unstable g-modes. 
The PMS and MS instability domain show a similar behaviour in 
the full Teff range, apart from a region near Pk ~ 6 days where 
we have young massive PMS models. To study this behaviour 
we focus on two non-adiabatic quantities 
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7500 7000 6500 

T 

^ eff 

Fig. 8. Theoretical IS for MS (light grey), MS with overshoot- 
ing (grey) and PMS (dark grey) j Dor stars compared to the 
observed y Dor IS (thick black lines - HS02). The thin dotted 
lines represent evolutionary tracks for models between 1 .4 and 
2.5M0. 



- the dimensionless work integral (W) defined as 



m — 



^3/2 /"/m(^r55)(r3-l)dm 
2 VGMw ^ (|^,|2 + t(t+\) |^,|2) dm 



(9) 



with R and M respectively the radius and the mass of the 
model, G the gravitational constant, w the frequency of the 
mode, s the specific entropy, the third adiabatic expo- 
nent, * defines the complex conjugate of a quantity, and 
and are respectively the radial and the horizontal com- 
ponents of the displacement. The work integral of a given 
mode gives the contribution of a region of the star to the ex- 
citation/damping of this mode. Regions where W increases 
contribute to drive the mode, while regions where W de- 
creases contribute to its damping. Finally, a mode is unstable 
if W > at the surface of the star. 

the eigenfunction dPjP that is linked to the amplitudes of the 
oscillations. Only regions where bPjP is large enough play 
a significant role in the driving/damping of the mode. 



CO 

3- 




7000 



Fig. 9. Period range of ^ = 1 unstable modes for all the computed 
PMS (dark grey), MS (light grey) and MS with overshooting 
(grey) models as a function of thek effective temperature. 



Before getting to the heart of this study, it is important 
to recall the excitation and damping processes that oc cur in 

} ' Dor stars. As developed in Guzik et al. (2000) and Dupre tet al.l 
20051) . the pulsations in y Dor stars come from a periodic mod- 
ulation of flux at the bottom of the convective envelope, there- 
fore the behaviour of the unstable modes is strongly linked 
to the depth of the CE. If two models have the same CE 
depth, they then have the same y Dor modes excitation capac- 
ity. Nevertheless, this is not sufficient to obtain the same pe- 
riod ranges of unstable modes. Indeed, another important pro- 
cess is the radiative darnping t hat occurs in the central layers 
of the star Dupre t et al.l ([2005') show that y Dor g-modes are 
unstable if and only if the excitation mechanism at the bottom 
of the convective envelope overcomes the damping mechanism. 
In the case of high-order g-modes, the eigenfunctions oscillate 
quickly in the ^-mode cavity, leading to high values of their sec- 
ond derivatives, which play a role in the work-integral expres- 
sion. This oscillation leads to a value of the radiative damping 
that is high enough to be more efficient than the y Dor driving 
mechanism. For low-order ^-modes, the absolute value of the 
eigenfunctions is significant in the ^-mode cavity compared to 
its value in the more superficial layers, which implies an efficient 
radiati ve damping in the ^-mo d e cavity. Within the asymp totic 
theorv.lVan Hoolst et al.](ll998h . lDz"iembowski etai](l200ll) . and 
iGodartet al.1 (12009') propose a simple expression of the radiative 
damping 77, which is proportional to 



lad - V V„rfA^gL 



dr. 



V Pr5 

We can rewrite ?; in a dimensionless form as 



^6 r-n 



-Ar ■ 



I{x)Ax. 



(10) 



(11) 



The radiative damping depends strongly on the matter stratifica- 
tion in the central layers of the star because of the factors 1 /r^ 
and A^. 



Illlllllllllllllllllllllllll III II III II 



0.5 1 1.5 2 2.5 

Period (days) 

Fig. 10. Period range for £ - I and £ = 2 unstable modes of PMS 
(black) and MS (grey) I.SMq models. 



5.1. Non-adiabatic properties of PMS and MS y Dor 

Fig. [To] presents the period domain of { = I and ( - 2 unsta- 
ble modes for two models of I.SMq. One is an MS model and 
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-0.2 




Fig. 11. Work integral (top panel) and eigenfunction (bottom 
panel) profiles for a PMS (black) and an MS (grey) modes with 



models explain the similarity of the period domain of unstable 
modes. 




Fig. 12. l(x) (oc radiative damping) versus the normalized radius 



the same dimensionless frequency. The vertical line represents for the 1 .8M0 PMS (full black line) and MS (dashed grey line) 



the bottom of the convective envelopes. 



models. 



Table 1. Model, angular degree, radial order, frequency and 
dimensionless frequency of the modes compared in the non- 
adiabatic section 





t 


k 


frequency 


dimensionless 








(c/d) 


frequency 


PMS I.8M0 


1 


-43 


0.566 


0.188 


MS I.8M0 


1 


-51 


0.576 


0.188 


PMS 2.1 Mo 


1 


-31 


0.662 


0.307 


MS 1.9Mo 


1 


-43 


0.637 


0.307 


PMS 2.1 Mo 


1 


-32 


0.637 


0.298 


MS 1.9Mo (a = 2.07) 


1 


-45 


0.614 


0.299 



the other a PMS one, both with the same effective temperature 
log Tetf = 3.85 (in the middle of y Dor IS, Fig. |2]- left panel). At 
this effective temperature, PMS and MS models present a similar 
instability domain. To understand the reasons for this similarity, 
we compare in Fig.[TT]the work integrals and the eigenfunctions 
corresponding to two unstable modes with the same dimension- 
less frequencjQ (see Table [TJ. It can be seen that for PMS and 
MS models the main driving mechanism leading to the excita- 
tion of the selected modes occurs at the bottom of the CE. Both 
models have the same CE depth (see Sect. O, hence their non- 
adiabatic functions present the same behaviour in the external 
layers. The only differences between PMS and MS eigenfunc- 
tions are located in the inner layers of the models and are due to 
their different central structures. Nevertheless, in both cases, the 
amplitude of ^ in the central layers is small, and the differences 
between both stellar structures have no impact on the excitation 
of the modes. 

Fig.[T2lillustrates the variation in the term I(x) (Eq. (fTTT i) for 
the I.8M0 PMS and MS models. In this case, I(x) (or what is 
the same, the radiative damping) behaves the same in both mod- 
els, excepted in their very internal layers. However, the impor- 
tant quantity to compare is the integral of I(x) on the radiative 
zone, i.e. the area under each curve. The close values of 77 in both 



' The oscillation frequencies (in c/d) in the MS and PMS 1.8 Mo 
models are close but not equal because the two models have slightly 
different radii. 



lllllflillMIIIIIIMIIilfPilTIIIIIIIIIIII 
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12 3 4 

Period (days) 

Fig. 13. Same figure as Fig. [10] for the 2.IM0 PMS and the 
I.9M0MS models. 



We focus now on the non-adiabatic behaviour of the 2.1Mq 
PMS and the 1.9Mq MS models, which are located close to the 
red border of the y Dor IS, where MS and PMS g-modes pe- 
riod ranges differ (see Fig. |9]l. Fig. [T3]presents the y Dor period 
ranges of these two models. While the number of unstable modes 
is the same for both models, the period domain corresponding to 
these modes is smaller for the MS model than for the PMS one. 
The different profiles lead to different periods and period spac- 
ing values (see Eq. ^ and Eq. As a consequence, the mode 
density is higher in the MS model than in the PMS one, and the 
instability range is narrower. 

As for the I.8M0 models (Fig. [TTT i. we compare in Fig. [14] 
the work integrals and eigenfunctions for two modes having 
the same dimensionless frequency. This mode is unstable in the 
PMS model but damped in the MS one. The two main differ- 
ences between these models are 

- the different driving zone depths due to different CE bottom 
locations, and 



7 



M.-P. Bouabid et al.: Theoretical seismic properties of PMS y Dor pulsators 




Fig, 14. Same figure as Fig. [TT] for the 2.1Mq PMS and the 
1.9MoMS models. 



- the important radiative damping that occurs in the inner lay- 
ers of the MS model due to a short wavelength oscillation of 
the eigenfunction in the g-mode cavity. 

To analyse the importance of each of these two factors on 
the damping of the MS mode, we displaced the CE bottom in 
the MS model by increasing the mixing length parameter from 
a = 2.00 to = 2.07. We obtain a I.9M0 MS model with the 
same effective temperature and the same temperature at the base 
of the CE as our 2.1Mq PMS model. Fig. [15] presents the same 
quantities and functions as Fig. [T3]for these two models. We still 
obtain an MS instability range narrower than the PMS one ow- 
ing the radiative damping which is much more efficient in the 
MS model (Fig. [16]- top panel). In Fig. [17] we present the I(x) 
profile of the two models discussed in Fig [14] We clearly see 
that the area under the MS curve is larger than the one under the 
PMS profile, especially in the inner regions where the MS model 
presents an important V^. The corresponding sharp peak in the 
profile leads to short wavelength oscillations of the eigenfunc- 
tions in the central layers, hence to a strong radiative damping of 
the modes. The decrease in the amplitude of the eigenfunction 
is so important that the driving that occurs at log T ~ 5.3 and at 
the bottom of the CE is not strong enough to drive the mode. 



6. Summary and conclusions 

In this paper we have presented the results of a theoretical study 
of the seismic properties of MS and PMS models inside the 
observational y Dor instability strip. In particular, we analysed 
whether the stellar structure differences are reflected in adiabatic 
and non-adiabatic seismic features, and whether it is possible to 
use seismic properties to distinguish between MS and PMS evo- 
lutionary phases of y Dor pulsators. 

The low-frequency g-modes that are excited in y Dor stars 
are mainly sensitive to the physical properties in the central re- 
gion of the star The more significant difference between PMS 
and MS models in the concerned mass domain is the presence of 
a chemical gradient in standard MS models, while PMS ones are 
almost homogeneous. The results from the comparison of adia- 
batic and non-adiabatic computations for both types of models 
are the following 




12 3 4 

Period (days) 

Fig. 15. Same figure as Fig. [10] for the 2AMq 
1 .9Mq MS models with the same temperature at 
their CE. 



PMS and the 
the bottom of 



0.001 




Fig. 16. Same figure as Fig. [TT] for the 2AMq 
1 .9Mq MS models with the same temperature at 
their CE. 



PMS and the 
the bottom of 



0.2 

X 



Fig. 17. Same figure as Fig. [H] for the 2.IM0 PMS and the 
1.9Mo MS models. 



- Due to the evolution of the convective core in PMS and MS 
phases, the value of the period spacing, resulting from the 
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first-order asymptotic approximation, changes as the stars 
evolve. Generally, (AP) is larger for PMS models than for 
MS ones, the maximum values being obtained when the 
convective core reaches its maximum size during the PMS 
phase. Nevertheless, close to the MS, PMS and MS models 
can have similar values of (AP); therefore it will not always 
be possible to distinguish between both phases of evolution 
on the base of only t he (AP > value. 

- As shown by Migli o et alj (2008), the //-gradient developed 
during the MS evolution leads to a periodic variation in the 
period spacing whose properties are related to the location 
and magnitude of that gradient. On the contrary, the almost 
homogeneous PMS models follow the first-order asymptotic 
approximation quite well and present a regular pattern of AP. 
Therefore, the regularity or variability of AP could be used 
as an indicator of the evolutionary state. 

- The PMS and MS instability strips overlap in the region of 
the HR diagram where models in both evolutionary states 
exist. 

- The period domain of unstable modes is quite similar for 
PMS and MS y Dor stars except close to the red border of 
the IS. In this region of the HRD, evolved MS models have a 
huge at the CC limit, which leads to a radiative damping 
that is much more important than in PMS models with quasi 
chemically homogeneous structures, and therefore the radia- 
tive damping aff'ects the oscillation properties of both kinds 
of models very differently. 

- The non-adiabatic computations also show a regular pattern 
of period spacing for high-order g-modes in PMS models, 
while a periodic dependence on the radial order is obtained 
for g-modes in MS models (see Fig.|7]l. 
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The comparison of these predictions with j Dor in young 
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